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Abstract 

We study spherically symmetric oscillations of electrons in plasma in the frame of classical 
electrodynamics. Firstly, we analyze the electromagnetic potentials for the system of radially 
oscillating charged particles. Secondly, we consider both free and forced spherically symmetric 
oscillations of electrons. Finally, we discuss the interaction between radially oscillating electrons 
through the exchange of ion acoustic waves. It is obtained that the effective potential of this 
interaction can be attractive and can transcend the Debye-Hiickel potential. We suggest that 
oscillating electrons can form bound states at the initial stages of the spherical plasma structure 
evolution. The possible applications of the obtained results for the theory of natural plasmoids 
are examined. 

PACS numbers: 52.35.Fp, 92.60.Pw, 74.20.Mn 

1 Introduction dilations of electrons in plasma [3]. This kind 

of plasma pulsation does not have a magnetic 

Theoretical explanation of the existence of self- field ( see Sec - EH below) and thus does not lose 

sustained spherically symmetric plasmoids is still energy for radiation. This fact can explain the 

a challenge in plasma physics [1]. A static bunch relative stability of spherical plasmoid generated 

of charged particles confined by its own inter- m natural conditions. 

nal forces does not seem to be stable. The hy- It iSj however, known that the frequency of free 

drodynamical pressure p and the magnetic pres- electrons oscillations in plasma cannot be less 

sure B 2 /8ir will try to expand a plasmoid (see t h an the Langmuir frequency. For an electrons 

Ref. [2J) making its existence impossible in the density in plasma of ~ 10 15 cm" 3 the Langmuir 

absence of external forces such as gravity, etc. frequency is about 100 GHz. It is rather difficult 

However a spherical plasmoid can be imple- to create a strong external field of such a high fre- 

mented in the form of spherically symmetric os- quency in natural conditions to generate a plas- 
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moid. Therefore one should point out a physi- 
cal mechanism that acts at the initial stages of 
the plasmoid evolution, during which one can en- 
counter a relatively low-frequency external field 
and which makes the plasmoid generation possi- 
ble. We suggest that plasma superconductivity 
can be one such mechanisms. 

The idea that a dense plasma can reveal su- 
perconducting properties was put forward in 
Ref. [1]. The superconductivity seems to play 
a key role at the stages of the plasmoid forma- 
tion. However it is difficult to find the physi- 
cal mechanism which would be responsible for 
the appearance of the superconducting phase in 
rather hot plasma. 

It is known that the formation of a bound 
state of two electrons, a Cooper pair, underlies 
the superconductivity phenomenon in metals. A 
Cooper pair is formed when the electrostatic re- 
pulsive field of an electron is shielded by the ef- 
fective attractive interaction due to the exchange 
of virtual phonons. Such a bound state is de- 
stroyed when the temperature of a metal exceeds 
a few kelvin degrees [5]. 

Although the temperature of plasmas in nat- 
ural conditions is far greater than a typical tem- 
perature of a superconducting metal one can find 
physical processes leading to the appearance of 
the effective attraction between electrons. A 
charged test particle, e.g. an electron, moving 
in plasma is known to emit ion acoustic waves. 
Therefore a test electron can be surrounded by 
a cloud of positively charged ions. Under cer- 
tain conditions this effective potential can screen 
the repulsive interaction between two electrons 
and result in the creation of a bound state (see 
Ref. [6]). This phenomenon, as well as the ex- 
change of dust acoustic waves, can lead to the 
effective attraction of dust particles in a dusty 
plasma [7j. Note that the formation of bound 



states of electrons in plasma due to the exchange 
of ion acoustic waves is analogous to the Cooper 
pairs formation [8j. 

In this work we study spherically symmetric 
plasma structures using the method of classical 
electrodynamics. Note that in Ref. [3] we con- 
sidered a spherical plasmoid as quantum oscilla- 
tions of electrons in plasma. Firstly, in Sec. 12.11 
we examine the electromagnetic potentials for 
the spherically symmetric motion of charged par- 
ticles and find a gauge in which the vector poten- 
tial is zero. Secondly, in Sec. 12.21 we study free 
and forced spherically symmetric oscillations of 
electrons on the basis of the system of equations 
of classical plasma hydrodynamics. In Sec. [3l us- 
ing the methods of Ref. [6] we calculate the scalar 
potential created by a test electron participat- 
ing in forced oscillations in plasma. We examine 
the conditions when the effective potential is at- 
tractive and consider the possibility of forming a 
bound state. Finally, in Sec. EJ we examine pos- 
sible applications of the obtained results to the 
description of natural spherical plasmoids. 

2 A description of spherically 
symmetric oscillations of 
electrons in plasma based on 
classical electrodynamics 

In this section we will study oscillations of elec- 
tions in plasma which have the spherical sym- 
metry. The method of classical electrodynamics 
will be used to describe this process. Firstly, 
in Sec. 12. 1| we will be interested in the various 
choices of electromagnetic potentials for such a 
system. Then, in Sec. 12.21 we will obtain the 
exact solution of the hydro dynamical equations 
which describes oscillations of electrons density 
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as well as the dispersion relation for these oscil- 
lations. 

2.1 Electromagnetic potentials in the 
system of radially oscillating par- 
ticles 

It is clear that the scalar potential of an electric 
field ip in the system making spherically sym- 
metric pulsations can depend only on the radial 
coordinate, <p(r,t). The vector potential A has 
only radial component which is also a function 
of only radial coordinate, A r (r, t). Due to the 
existence of the spatial dispersion in plasma, os- 
cillations of charged particles should vanish on 
big distances from the center of the system, i.e. 
A r (r,t) — > at r — > oo. 

Suppose that we have found the potentials ip 
and A. Now we can make the gauge transfor- 
mation, A' = A + V/ and tp 1 = ip-\- (l/c)df/dt, 
where f(r,t) = f°° A r (r',t)dr'. This gauge 
transformation does not change the electric field. 
Note that the magnetic field is identically equal 
to zero for a spherically symmetric motion of 
charged particles, B = V x A = 0. For the 
chosen function / we obtain that the vector po- 
tential can be eliminated in all the space. Now 
the electric field has the form, E = — Vy>'. In the 
following we will omit the prime in the definition 
of the scalar potential. 

2.2 Classical plasma hydrodynamics 
description of a spherical plas- 
moid 

In the first approximation we suggest that only 
electrons participate in oscillations of plasma 
since the mobility of ions is low. In the absence of 
collisions and other forms of dissipation the sys- 
tem of the hydrodynamic equations for the de- 



scription of plasma oscillations can be presented 
in the following way (see Ref. [9]): 



dn e _ . . 



-Vp, 



dv , „ s e „ 

+ v-V v = E 

at m 

V • E = -Aire(n e - n») + 47rp ext (r, t), 

dB 

— = 4vren e v, 
at 
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where n e is the electrons density, rij is the ions 
density, p is the plasma pressure, m is the mass 
of the electron and e > is the proton charge. In 
Eq. ([T]) we include the possible external source 
p ex t and take into account that the magnetic field 
is equal to zero (see Sec. 12. ip . 

Supposing that nj = no and the deviation of 
the electrons density from the equilibrium value 
is small, n e — uq = n <C no, we can linearize 
the system (TTJ) and obtain the single differen- 
tial equation for the perturbation of the electrons 
density, 

d 2 n 2 1 / dp s 
dt 2 e m \ dn i n in 



2 4vren 
V n = p^, (2) 



where uj e = \J ' A / Ke 2 no/m is the plasma frequency 
for electrons and (dp/dn)o is the derivative taken 
at n e = no- The latter quantity depends on the 
equation of state of electrons in plasma. 

In the absence of the external source the spher- 
ically symmetric solution of Eq. ([2]) has the form 



n(r, t) = Acos(ut) 



sin 7r 



(3) 



where A is the constant chosen to satisfy the 
condition |n| no- The frequency of oscillations 
uj and the length scale parameter 7 are related 
by the following identity: 

1 (dp 
m V dn 



wf + 



7 



(4) 
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which shows that free oscillations with u > oj e 
exist in plasma. However, if p cxt ~ cos £lt, it is 
clear that forced oscillations with Q < u e can be 
also excited. 

In Ref. [3] we studied spherically symmet- 
ric oscillations of electrons in plasma using the 
quantum mechanical approach. In that paper, 
we solved the non-linear Schrodinger equation 
for the wave function normalized on the num- 
ber density of electrons, |^(r, i)| 2 = n e (r,t). It 
was obtained that in the spherically symmet- 
ric case the density of electrons has the form, 
n e (r, t) = hq + Acos(cot) sin(7r)/r + • • • , i.e. is 
similar to Eq. ([3]). However the dispersion rela- 
tion in Ref. [3] was different from Eq. ([3]). More- 
over quantum oscillations of electrons reveal the 
typical size of the system, where the most inten- 
sive oscillations happen, L = tt/j = TTy / h/2mu} e , 
at the critical frequency cj = 2u e . On the con- 
trary, if we use the classical electrodynamics 
method, the parameter tt/j, in principle, can be 
arbitrary. Of course, Eq. Q is valid only for 
quite long waves, when 7 <C k e (see Ref. [10]), 
where k e is the Debye wave number (see the def- 
inition in Sec. [3]). 



3 Formation of bound states of 
electrons at the initial stages 
of the spherical plasmoid evo- 
lution 

Any particle in plasma has a two-fold life. On 
the one hand it is a test particle moving through 
plasma and interacting with the whole plasma 
rather than with separate plasma particles. On 
the other hand any test particle is a part of 
plasma and hence it contributes to self-consistent 
electromagnetic fields in plasma. In this sec- 



tion we will use the method of test particles (see 
Refs. [6, 8j) to calculate the effective potential of 
a radially oscillating electron. 

Let us study the electric field created by elec- 
trons participating in spherically symmetric mo- 
tion considering each electron as a test particle 
with the charge q. Each of the test particles is 
taken to interact with the rest of hot electrons, 
having the temperature T, and with cold ions. 
The permittivity of this plasma has the form, 

^)=i+^Y-f-^y, (5) 



k 



u + iO 



where k e = \J Aim^e 2 /T is the Debye wave num- 
ber and Ui = \J Air(Ze) 2 no/M is the plasma fre- 
quency for ions with the mass M and the charge 
Ze, Z is the degree of the ionization of an ion. 

We suppose that a test particle makes har- 
monic oscillations around the point ro with the 
frequency Q and the amplitude a: r'(t) = tq + 
a sin tit. To find the electric potential created by 
one of the charged particles we should account 
for the Maxwell equation for the electric dis- 
placement field, V • D = 47rq<5 3 (r — r'). Express- 
ing the electric field as E(k, u) = — ik</?(k, cj), we 
obtain the scalar potential of the system in the 
form 

ip{k,u)) 



k 2 e{k, 



We recall that the vector potential can be taken 
to be equal to zero for the spherically symmetric 
system (see Sec. [HQ. 

Let us decompose the exponential factor in the 
integrand of Eq. © using the series of Bessel 
functions of the n-th order, 



-ik-asin£M 



ale 



(7) 
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Finally on the basis of Eqs. © and ([7]) we obtain Taking into account the value of the following 
the scalar potential in the form integral: 



+00 



"2vr 2 



^ ! do;d 3 ke- iwt+ik -( r - r °) 
=-00 J 

J n (k ■ a) 



x 5(ui — nCl) 



fc 2 e(k,w)' 



(8) 



To analyze Eq. (Jf)]) we present the reciprocal 
of the permittivity © in the following way: 



k 2 



e(k, oj) k 2 + k 2 \ oj 2 — oj 2 



(9) 



where u a = kcoi j \Jk 2 + k 2 is the dispersion re- 
lation for ion acoustic waves. 

Using Eq. ([9|) we present ip in Eq. ([8]) as a sum 
of two terms, (p = pr> + <pw-> where 

+00 „ 

^ D ( r)t ) = _i_ ^ / d 3 ke -inm+ik.(r-r ) 



n=— 00 

Jn{k- a) 

k 2 + k 2 ' 



(10) 



is the analog of the Debye-Hiickel screening po- 
tential and 



ipw 



+00 . 

^=^2 E y d3ke 



— inf2t+ik'(r— ro) 



n=— 00 

(nft) 2 - lo 2 k 2 + k 2 ' 



(11) 



is the wake potential due to the emission of ion 
acoustic waves 

To study the behaviour of the potentials pr> 
and ipw we choose the specific coordinate system 
with ro = and a = ae z , where e z is the unit 
vector along the z-axis. We also decompose the 
vectors k and r using the cylindrical coordinates: 
k = (k p ,k z ,(f>k) and r = (p,z,(f>). 



d^e ifc " cos ^-^ = 27rJo(M> (I 2 ) 
we rewrite the potential ^ in the form, 

<p D {p,z,t) =2q V(-l) n / k p dk p 



n=0 

e 



kj + k 2 



x |/ 2n (a^Jkj + kfj cos[2nfti] 

± ^2n+l { a \J^ 2 p + fcf) 

[(2n+l)J2t]}, (13) 



x sin 



where I n (x) = (i) _n J ra (ix) is the Bessel function 
of the imaginary argument. In Eq. (|13p the '+' 
sign stands for z > and ' — ' for 2 < 0. 

The electromagnetic field of a charged linear 
oscillator in vacuum, e = 1, was studied in 
Ref . |11| . The scalar potential found in that book 
is different from Eq. (|13p since in Ref. [11] the 
problem of radiation of a linear oscillator was 
considered and the Lorentz gauge for potentials 
was used. We study the case of the spherically 
symmetric motion of plasma. As we demon- 
strated in Sec. 12.11 such a system does not have 
any magnetic field and thus cannot emit radia- 
tion. In our situation it is more convenient to 
use the gauge in which A = 0. That is why we 
get a different expression for the scalar potential. 

It is interesting to analyze Eq. (|13p in the 
static limit. For the test particle at rest, i.e. 
at a — >■ 0, only the term with the function Iq(x) 
survives. It is possible to show that in this limit 
Eq. (fTBj) transforms to ip£)(p,z,t) = (q/r)e~ ker , 
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where r = \J ' z 2 + p 2 . We can see that one recov- 
ers the usual form of the Debye-Hiickel potential. 
This analysis justifies our definition of </?£>. 

Now we study the wake potential (fTT|) using 
the same technique as for iff). We present the 
value of ifw in the cylindrical coordinates as 



-2q 



kpdkp 



U2 _i_ h.2 



x I (a^Jkf+k^j 
+ 2q 



oo 
n=l 



1 

(-1)" 
1 



0J7 



007 



nVL) 2 
dxJo(k n p\/ 1 — x 2 ) 



ke + JO 

x J n (ak n x) sui{k n \z\x — nflt), (14) 



where 



k n — k e 



nfl 



007 



(15) 



The upper multiplier in Eq. (jlij) corresponds to 
z > and the lower one - to z < 0. Comparing 
Eq. (I13p and Eq. (|14h we can see that the terms 
containing the integrals of Iq(x) cancel in the 
sum of the two potentials. 

To derive Eq. (fT4"J) we suggest that nfl < cjj, 
i.e. the wake potential appears only for non- 
rapid oscillations. For example, we can consider 
the forced oscillations of electrons in plasma de- 
scribed in Sec. l2.2l and suppose that Q is a bit less 
than LOi. It means that only the first harmonic 
is excited. 

Let us study Eq. (JT2J) at the line of the test 
particle oscillations, p = 0, and at big distances 
from the test particle, \z\ 3> a. Putting n = 1 
we obtain from Eq. (114H 



<pw(z,t) ~ T<7- 



tt 2 



oot 



n 2 



cos(ki\z\—Clt), (16) 



where the '— ' sign stands for z > and '+' for 
z < 0. 

Returning to the general Eq. (fl~3j) for </?£> we 
also rewrite it for /j = as 

2q ' 

\z 



^(z,f)=^(-l)» / dxe~ x 



e-^^m Jo (k pP ) 



n=0 



|-^2n f cos[2nf2t] 



2n+l 



-.T 



x sin[(2n + l)fit]j. (17) 

Comparing Eqs. (|16|) and (|17p we can see that 
at large distances the non-Coulomb wake poten- 
tial transcends the Debye-Hiickel potential. For 
example, the term with the function I\{x) in 
Eq. (fTT|) at \z\ 3> a and k e a <C 1 has the form 



T ^ e -M*l sin(m ), 



(18) 



which has much smaller value than the wake po- 
tential ()16p at \z\ > l/k e . The terms which con- 
tain the functions I n (x) with n > 1 will give con- 
tributions to ipr> smaller than that in Eq. (JTHJ) . 

The wake potential is attractive when 
cos(/ci|z|— fit) > for z > and when cos(/ci|z| — 
fit) < for z < 0. We recall that we study 
the radial pulsation of plasma. Thus the coordi- 
nate z coincides with the radial direction. Sup- 
pose that a test electron attracts another elec- 
tron which is situated at the distance d away the 
center of the system during the certain period of 
time, Aii = ir/fl. During the next half a period 
of the wake potential variation, At2 = n/fl, the 
same test electron will attract a different elec- 
tron situated at the same distance d but closer 
to center of the system. It means that a test elec- 
tron can always attract some charged particles in 
plasma. 
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To form a bound state with another electron 
in plasma the energy of interaction of a test elec- 
tron eipw should be greater than the total energy 
of its oscillations, E osc = ma 2 Q 2 /2. Suppose 
that two electrons are at the distance d = l/k e 
and the amplitude of oscillations a = O.ld <C 
d. Studying the plasma consisting of electrons, 
with temperature T = 10 3 K and number den- 
sity no = 10 15 cm -3 , as well as of singly ion- 
ized nitrogen atoms we obtain that the ratio 
| eipw | /E osc > 1 if fi > lCP 4 u>j. Taking into 
account that Wj ~ 10 10 s^ 1 we obtain that the 
frequency of the forced oscillations should be in 
the region 10 6 s _1 < < 10 10 s -1 . Thus we get 
that the effective attraction can take place in the 
atmospheric plasma for the reasonable frequen- 
cies of an external field. 

4 Discussion 

We studied electrons oscillations in plasma us- 
ing classical electrodynamics. It is found that 
spherically symmetric oscillations are possible in 
the classical case, although the dispersion rela- 
tion ([3]) is different from the previously found 
for quantum oscillations [3]. We suggest that 
the radial pulsation of plasma underlies the rare 
atmospheric electricity phenomenon called a ball 
lightning (BL) [I]. 

Radial oscillations of electrons in plasma de- 
scribed in the present work are unbounded and 
occure in all the space. It is the case when a plas- 
moid, BL, propagates far away from any external 
surface. The situation changes when one consid- 
ers spherically symmetric oscillations of plasma 
in a cavity inside a dielectric material. This kind 
of oscillations was studied in Ref. [12] in pres- 
ence of external electric and magnetic fields. The 
exact system of non-linear equations describing 



oscillating spatial patterns was obtained and an- 
alyzed both analytically and numerically. The 
solutions obtained in Ref. [12] can correspond to 
BL passing through microscopic cracks in a di- 
electric material, e.g. glass. There are several 
reports of such events collected in Ref. [13] , The 
problem of the interaction of BL with external 
materials should be analyzed in out future works. 

As we demonstrated, in order to generate 
spherically symmetric oscillations, one should 
excite them with rather high frequencies, Q > u e 
in the classical case (or £1 > 2u e in the quan- 
tum case). The electron plasma frequency is 
(10 12 - 10 13 )s~ 1 for n = (10 15 - 10 17 )cm~ 3 . 
Such high frequencies are very difficult to obtain 
in any natural conditions. 

We showed in Sec. 12.21 that forced oscillations 
are possible with frequencies less than cj e . Of 
course, forced oscillations will decay as soon as 
the external force is switched off. There should 
be a mechanism which provides the smooth tran- 
sition from the generation regime of a spherical 
plasmoid with the external harmonic source hav- 
ing Q < co e to a regime with self-sustained oscil- 
lations having Q ~ uj e . 

We suggest that this mechanism could be a 
formation of bound states of electrons in plasma. 
The motion of bound states of electrons can re- 
sult in the appearance of a superconducting state 
of plasma inside a spherical plasmoid. Previ- 
ously the idea that the superconductivity can ex- 
ist in plasma was put forward in Ref. [13] • With- 
out the existence of superconducting phase elec- 
trons participating in radial oscillations will lose 
their energy very quickly, because of the various 
friction mechanisms, and will recombine into the 
initial neutral gas. 

In Sec. [3] we obtained that a test electron os- 
cillating in plasma with the frequency cjj 
would emit ion acoustic waves. Thus the test 
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electron appears to be surrounded by a cloud 
of 'phonons' that shield its repulsive potential. 
Under some conditions the effective interaction 
between the test electron and other electrons in 
plasma turns out to be attractive. This process 
can lead to the formation of bound states of elec- 
trons. 

Note that for the first time the role of 
'phonons', or acoustic waves, for the descrip- 
tion of the stability of a spherical plasmoid was 
discussed in Ref. [15] . However, in that work, 
the phonons exchange between ions and neutral 
atoms was considered in frames of the quantum 
theory. The applications to the obtained results 
to the theory of BL were also studied in Ref. [15] . 

The existence of the superconducting phase 
inside BL was previously proposed in Ref. [JJ. 
However in these papers the dense plasma of 
BL was already supposed to be in the supercon- 
ducting phase and the phenomenological conse- 
quences of this phenomenon were described. No 
physical mechanisms for the formation of the 
superconducting state were proposed. In the 
present work we suggest that the exchange of 
ion acoustic waves created in spherically sym- 
metric oscillations of electrons in plasma results 
in the effective attractive potential between elec- 
trons. This mechanism can lead to the formation 
of bound states of electrons and possibly to a su- 
perconducting state inside a spherical plasmoid. 

After one switches off the low frequency ex- 
ternal field, bound states of electrons will be de- 
stroyed. We assume that during the supercon- 
ducting stage of the plasmoid evolution a source 
of internal energy should appear. It was sug- 
gested in Ref. [16] that nuclear reactions can 
serve as the energy source of BL. The recombina- 
tion of charged particles will be compensated by 
the processes of ionization and creation of new 
charged particles owing to this internal source of 



energy of BL. 
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